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Abstract 

Renormalization group transformations as discussed recently in deriv- 
ing fixed point actions are used to analyse the vacuum structure near to 
the deconfinement temperature. Monte Carlo configurations are gener- 
ated using the fixed point action. We compare equilibrium configurations 
with configurations obtained by inverse blocking from a coarser lattice. 
The absence of short range vacuum fluctuations in the latter does not in- 
fluence the string tension. For the inversely blocked configurations we find 
the following: (i) the topological susceptibility Xtop is consistent with the 
phenomenological value in the confinement phase, {ii) xtop drops across 
the deconfinement transition, (iii) density and size of instantons are es- 
timated, (iv) the topological density is found to be correlated to Abelian 
monopole currents and (v) the density of spacelike monopole currents be- 
comes a confinement order parameter. 
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1 Introduction 



Not much is known from first principles about the topological structure at the 
thermal phase transition in Quantum Chromodynamics. How this structure 
changes at the deconfincment phase transition is not even clear in pure gauge 
theory. What we mean by "topological structure" is the characterization of 
Euclidean field trajectories in terms of instantons [|l|, || and Abelian monopoles 
[|[ ^ and their interrelation. Changes of this structure are expected to be the 
common mechanism for deconfincment and restoration of chiral symmetry. Un- 
derstanding this mechanism is important for understanding the QCD vacuum. 

There was a wave of attention recently for improving lattice actions for gauge 
fields 1^ (and fermions). Part of the interest was motivated to minimize finite 
lattice spacing corrections in simulation results obtained on coarse lattices. One 
particular approach towards improvement was guided by the desire to construct 
a perfect action being a fixed point action |^ with respect to some real 
space renormalization group (RG) transformations. In this context the concept 
of inverse blocking came up which is the essence of a new method, recently 
proposed by DeGrand et al. ^. In particular it improves the determination 
of the topological charge of a lattice configuration. 

In the present paper we want to explore the capability of this method, to 
resolve more details and different aspects of topological structure. In the neigh- 
bourhood of the deconfincment transition of the pure SU (2) gauge theory, we 
want to study the space-time distribution of topological charge, correlations of 
the density with itself and with Abelian monopole currents. 

The topological charge density, when defined for quantum fields, is obscured 
by ultraviolet fluctuations which must be removed in some way. Cooling [ p^ has 
been invented for this purpose long ago. According to the Abelian dominance 
scenario at large distances El Abelian magnetic monopoles are closely related 
to confinement. The density and anisotropy of these magnetic currents 
are studied in the maximally Abelian gauge fl^ . These currents are quantized 
by construction This does not exclude the possibility that some part of 

these currents are lattice artefacts which do not correspond to long distance 
physics. Correlations indicating a close relation between topological charge and 
Abelian magnetic monopole currents have been studied, not addressing this 
question. They have been studied so far for particular instanton configurations 
and, as far as the quantized vacuum is concerned, the cooling method 
has been employed [p^ . However, cooling is known to change the topological 
structure and has to be used with care. 

A closer examination shows that the method of Refs. ||] gives not only 
the topological charge but also permits a reasonable definition of the topolog- 
ical density of generic, equilibrium lattice fields with minimal effort. A better 
resolution of topological and monopole structure is possible without recourse 
to cooling. While the method is applied to extract the topological structure, 
the monopole activity is changed, too, roughly in proportion to the topological 
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activity (to be defined below) . We shall see that part of it becomes an order pa- 
rameter of the confinement /deconfinement transition, according to the alleged 
role monopoles play for confinement. 

According to the scope of our present study, let us sketch the background 
of this approach |^ restricting the discussion to the case of pure gauge field 
theories. The starting point is a block spin transformation (here of scale factor 
two) which maps a fine lattice configuration of link matrices Ux^fj, to a coarse 
lattice configuration with link matrices (each single link being an element 
of the gauge group SU{N)). The second ingredient is a general type of action 
Spp, written on both lattices in terms of Wilson loops (evaluated in various 
representations). The classical perfectness of this action is guaranteed if there 
exists an inverse blocking transformation V U with respect to which the 
action saturates the inequality 

SFpiU) + ^iT{U,V)>SFpiV) . (1) 

Here, T{U, V) is some non-negative functional related to the blockspin trans- 
formation and K is a kind of Lagrange multiplier. Once the action and other 
parameters have been determined, this inequality leads to the target configura- 
tion U of the inverse blocking by minimization of the 1. h. s., the fixed point 
action with respect to U, with a constraint derived from the coarse configu- 
ration V. Thus, U is found as the result of an iterative relaxation process 
and cannot be explicitely expressed. We call the mapping one step 

of inverse blocking. Similar to blocking, inverse blocking can be iterated. This 
leads to a sequence of interpolating configurations on finer and finer lattices. 

Without giving details, we mention that the saturated inequality (|^) is con- 
structive for finding the perfect action itself. It has been used Q to optimize the 
ratios between the couplings in front of various loops in different representations 
within the fixed point action S'fp, as well as the other constants related to the 
block spin transformation. 

After this machinery has been set up, the procedure of inverse blocking re- 
mains a useful tool of analysis. Doing a sequence of inverse blocking steps one 
would be able to relate a continuum configuration to any given lattice configu- 
ration. The topological charge of this continuum configuration is the topological 
charge related to the given lattice configuration (the classically perfect topologi- 
cal charge) . This prescription is an important conceptual step forward compared 
with previous techniques of measuring topological charge on the lattice, even- 
tually used in combination with cooling. 

Inverse blocking has been used before ||, only to measure the topological 
charge of lattice configurations, by evaluating the Phillips-Stone (PS) charge 
[ p!?! just on the next finer level. The success indicates that one step of inverse 
blocking can be sufficient to have access to the classically perfect topological 
charge corresponding to the continuum by means of a lattice definition. Other 
lattice algorithms besides PS could also have been used, for instance in order 
to test the uniqueness of charge. We will use in this paper a geometric charge 
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due to Liischer [I^, ^ and the old plaquette based (field theoretic or naive) 
definition |2^, |21 1 . Both have the advantage to localize the topological density 
in space-time which makes them suitable for our present purpose. 

The first concern in Ref. |8| has been that action and topological charge of 
artificially constructed lattice instantons should be independent of size, if the 
perfect action is used and the topological charge is measured as described. This 
was demonstrated for instanton configurations with p > 0.8 • a. There were 
no dislocations, with action S/Sj < 6/11 (cntropic bound, instanton action 
Si = 27r^) and non-zero topological charge, among them. However, the perfect 
action and the new definition of topological charge are useful also at finite /3. In 
the second paper ||] a (truncated) fixed point action has been successfully tested 
for scaling of the deconfinement temperature, of the torelon mass, of the string 
tension and of the topological susceptibility. For the latter, it has been stressed 
that inverse blocking has been indispensable to achieve reasonable scaling. This 
has nourished the hope that the fixed point action and the inverse blocking 
prescription for the topological charge can avoid dislocations which otherwise 
(for instance in the case of Wilson's action) are known to spoil the continuum 
limit of the topological susceptibility |Q . 

Even the temperature dependence of the topological susceptibility has re- 
mained controversial until recently. Armed with the new method, the topo- 
logical susceptibility has been measured across the deconfinement transition |^ 
without the ambiguities of the cooling method ||2^, . Inverse blocking is not 
the only technique to avoid cooling. We only want to mention the use of an im- 
proved lattice operator of the topological density |Q . This approach has given 
a much more discontinuous drop [ p6| at the deconfinement transition than the 
inverse blocking method |^ . 

Our paper is organized as follows. In section 2 we give details on the fixed 
point action, on blocking and inverse blocking. In section 3 we present results 
on the string tension. For symmetric lattices this serves to calibrate the lat- 
tice step size a(/3) in accordance to simulations using the fixed point action. 
The finite temperature string tension is computed (in the confined phase) from 
Polyakov line correlators. These results are indicating that the string tension 
is largely conserved by inverse blocking. Section 4 contains our results on the 
gross topological features. The topological susceptibility obtained on inversely 
blocked configurations in the confined phase is quantitatively reasonable. The 
temperature dependence of the topological susceptibility across the phase tran- 
sition is studied. Analysing inversely blocked configurations we observe, apart 
from the decay of the topological susceptibility, the complete suppression of spa- 
tial Abelian monopole currents in the deconfined phase. In section 5 we present 
results on topological charge density-density correlations and on monopole- 
instanton correlations. We visualize clusters of topological charge together with 
their accompanying monopole currents for typical confining configurations. In 
closing, directions for further work are pointed out in section 6. 



4 



2 Fixed Point Action and Inverse Blocking 
2.1 Action and Blocking 

The simplified fixed point action Q is parametrized in terms of only two types 
of Wilson loops, plaquettes Uc^ = U^^^^u (type Ci) and tilted 3-dimensional 
6-link loops (type C2) of the form 

It contains several powers of the corresponding linear action terms as follows 
Sfp{U) = J2 EE^(*'J')(1 - ^ UC'V ■ (3) 

type i d j—1 

The parameter of this action have been optimized in Ref. |^ and are reproduced 
in Table 0. 



w{i,j) 


J = l 


J = 2 


J =3 


J = 4 


i = 1 (plaquettes) 
i = 2 (6-link loops) 


.3333 
.08333 


.00402 
.0156 


.00674 
.0149 


.0152 
-.0035 



Table 1: Weight coefRcients of the simphGed fixed point action 

The notion of blocking is recalled here in order to define some notations. 
Fat links I4_^ are defined for Xi, = odd for all v = 1, 2, 3, 4 only. These are the 
nodes of the blocked lattice.^] The blocking transformation expresses a fat link 
in terms of the product of the two links along the fat link plus a sum over six 
rectangular 4-link staples passing by the fat link: 

V~' block jj Tj 



+ E] ' {yx,uUx+v,^iUx+v+ii.4iUx+2ii,v ('^) 



The actual fat link Vx^^ is then obtained by normalization to SU{2): 

Vx,, = (5) 
det (Vx^fj, 



^The blocked lattice could be defined as well shifted by one fine lattice step along one or 
more directions. 
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The following blocking parameters have been recommended for this type of RG 
transformation in Ref. |Q 

^block ^ 

^uock ^ i-6c^'°"'= . (6) 



2.2 Inverse Blocking 

Inverse blocking consists in searching the constrained minimum with respect to 
the fine links U of an extended action 

S,,t{U)^SFp{U) + KT{U,V) , (7) 

which includes, besides of the fixed point action Spp, the blocking kernel 




with t\ clS Oil Lagrange multiplier. For better readability, we have dropped the 
labels x,iJi of the fat links 14. ^ (describing the blocked configuration) and of the 
Vx.fj.- The latter are expressed in terms of the fine J7's exactly as in (jj) and are 
not related to Vx.^i in the present context. 

For comparison, cooling the gauge field configuration U means to minimize 
the new action with k = 0. The result of cooling is no more influenced by 
the coarse configuration V in the process of relaxation. The unconstrained 
minimum [/ is a classical solution on the lattice with respect to Spp. It can 
only depend on the initial configuration if that determines a non-trivial basin 
of attraction. 

Concerning the organization of measurements, we differ from DeGrand et 
al. They have simulated on the coarse lattice, using the fine lattice only for 
measuring the topological charge on the inversely blocked configurations. The 
coarse lattice was their reference scale. We carried out simulations on the fine 
lattice (at some lattice scale a) chosing a /3-value to produce a Monte Carlo 
ensemble of gauge field configurations U^'~' according to the fixed point action. 
Most of these configurations have been used only for comparison as described 
below. Another part of them (in a separate run with the same [3) has been 
used to produce an ensemble of coarse lattice configurations V by blocking the 
fine lattice configurations. Due to variance reduction it was possible to restrict 
this second ensemble to considerably less configurations than the other. At 
scale 2 • o (and bigger) the coarse ensemble is expected to describe the same 
physics as the fine ensemble. We could have asked to which (3' the coarse lattice 
configurations correspond. However, for the sake of clarity in the comparison 
we always refer to the original /3-value and a as our reference scale. The coarse 
lattice configurations V have been analyzed with respect to their topological 
properties by inverse blocking, V U^^ . 
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To create the coarse configurations by blocking from another equilibrium en- 
semble is no matter of principle. It is, however, an advantage to keep the Monte 
Carlo configuration on the fine lattice in order to use it as a start configuration 
for the constrained minimization. Due to this particular circumstance the in- 
versely blocked configuration U^^ can be imagined as a smoothed copy of the 
original configurations U^^'~^ . Nevertheless, it can exhibit no other topological 
structure than present in the blocked configuration V (which influences the 
relaxation through the blocking functional kT(U, V)). In what follows we shall 
denote the combined application of one step of blocking followed by one step of 
inverse blocking as smoothing. 

Classical configurations are particular. Instanton solutions on the fine lattice 
are always reproduced under smoothing. Their topological structure is the same 
on the coarse and the fine lattices. Not only the action is recovered but also 
the shape of the topological charge distribution and the location of monopole 
currents. We have produced a set of almost classical configurations by cooling 
from hot starts using the simplified fixed point action (|^). This set has then 
been used for testing our code of the smoothing algorithm. If one blocks a 
classical solution, one finds Sfp{U) = Sfp{V). 

However, this technique is intended to be applied to non-classical configu- 
rations. Before entering our investigations, we have tested, creating fine lattice 
configurations U^'~' (Monte Carlo generated with Spp at various /3-values) and 
blocking them into coarse configurations V, which K-value must be used in the 
extended action (|^) driving U from J/*^*^ to U^^ . The relaxed configurations 
U^^ should turn the inequality (|l|) into an equality relating U^^ to V. Un- 
der these circumstances we were not able to confirm the Lagrange multiplier 
K = 12 recommended in Ref. (for the given type of RG transformation in 
SU{3)). Instead, we have found that a choice k = 5.15 (slightly depending on 
P) is suitable, throughout the /3 range under study, to guarantee the equality 
in ([l]) within an accuracy of one to three per cent. We have adopted the fixed 
point action but decided to use our empirical value for k. In principle, k and the 
parameters of the fixed point action are closely tied together. They should be 
optimized w. r. t. (Q) using an ensemble of configurations created with exactly 
this action. This is a much more ambitious project than it has been pursued 
until now. The authors of Ref. are continuing this program.^ This whole 
discussion puts a question mark on the exact relation between the action and k. 
Also in view of this, we consider the results of the present paper as exploratory. 

Knowing the extended action (|^) we can distinguish three types of links 
on the fine lattice. Links along the fat links and links lying inside coarse pla- 
quettes contribute to V and are updated (relaxed) with respect to Sext- All 
other fine links are not subject to constraints and are updated with respect 
to Spp alone. A typical relaxation process is shown in Fig. 0. The upper 
curve demonstrates how the extended action Sfp{U) + n T{U, V) decays versus 

^private communication 
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Figure 1: Relaxation towards the inversely blocked configuration vs. cooling for 
a typical Monte Carlo configuration generated at P = 1.5 on a 12^ x 4 lattice 



number of iterations. The upper horizontal line marks the action of the original 
Monte Carlo generated configuration Sfp{U^'~'), the lower horizontal the action 
of the blocked configuration Sfp{V). After V has been obtained from [/*^*-^ 
by blocking and one has plugged in U'^^'~'' into the extended action Sext{U), 
T{U, V) vanishes by construction. Therefore the minimization of Sext starts 
at the level of Sfp{U'^^'~^). After the constrained minimum configuration U^^ 
has been found, the extended action is shared between the fixed point action 
Sfp{U^^) and the functional n T{U^^,V). In Fig. |l| the history of Sfp{U) 
is separately recorded by the curve in the middle. The lowest curve shows how 
Sfp{U) would evolve under cooling (with k put to zero). After 20 iterations the 
unconstrained relaxation runs into configurations U with a fixed point action 
Sfp{U) one order of magnitude below Sfp{U^^) and the action decays further. 
On the other hand, this demonstrates that inverse blocking does not need huge 
numbers of iterations to achieve convergence if a good start configuration is 
known. A particular relaxation scheduling, for instance a modification of the 
action during the first iterations, is not necessary. 
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3 String Tension at T = and Finite 
Temperature 



The inverse coupling constant (3 enters through the Gibbs weight 

Z = J DUexp{-(3SFp{U)) (9) 

into the creation of fine lattice configurations. It defines the scale through the 
lattice spacing a{j3). We wanted to compare the naive charge with Liischer's 
charge which is, similar to the fixed point action, very CPU intensive. There- 
fore we have restricted ourselves to a fine lattice of size 12'^ x 4 for the finite 
temperature part of these explorative investigations. 

For i( — 4 the deconfinement transition was known to occur at f3c = 
1.575(10) §. This has been established by the crossing of the Binder cumulants 
of the Polyakov line for several spatial lattice sizes. We have not attempted here 
to improve the localization of the deconfining transition or to do a finite size 
analysis. Rather than doing that, we have simulated at /3-values safely lower 
(/? = 1.40, /3 = 1.50 and /3 = 1.54) and higher (/3 = 1.61 and /3 = 1.80) than the 
deconfinement transition. 

3.1 Lattice Spacing from the Zero Temperature 
String Tension 

In order to measure the lattice spacing a at each of the /3-values to be used in 
the finite temperature simulations we have complementarily measured the zero 
temperature string tension on a 12^ lattice. No topological features have been 
analysed in this part of our calculations. The string tension has been obtained 





1.40 


1.50 


1.54 


1.60 


1.80 


Monte Carlo 
inv. blocked 


0.445(8) 
0.228(4) 


0.245(3) 
0.182(1) 


0.193(3) 
0.156(3) 


0.124(2) 
0.108(3) 


0.0380(9) 
0.0299(1) 


Monte Carlo 
inv. blocked 


152(3) 
77.9(10) 


135(2) 
100.5(10) 


129(2) 
104.5(20) 


111(2) 
96.8(20) 


91.7(19) 
71.2(20) 



Table 2: String tension a (above) and (t/A| (below) at zero temperature from 
fuzzy Wilson loops in axial gauge on a 12* lattice (with two-loop expression for 
aiP)AL) 

for Monte Carlo configurations in the axial gauge from smeared Wilson loops. 
The same has been measured on smoothed configurations for comparison. The 
axial gauge is enforced in chronological order time slice by time slice (except for 
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-2 
smear 



the last one) after a certain direction has been chosen as Euclidean time. This 
is done putting each temporal link U4{x) equal to unity by applying a suitable 
gauge transformation g{x + 4) at the end of that link. Smearing |^ has been 
applied to configurations in this gauge. This is an iterative operation (without 
change of scale) replacing each spatial link U by through a procedure 

similar to (jj) (since without change of scale, with only one step in /x = 1,2,3 
direction and with only spatial staples involved, z/ = 1, 2, 3 v fi). In the axial 
gauge, timelike parts of the Wilson loops are mostly equal to unity (if they 
do not contain links from the last to the first time slice). Smearing effectively 
replaces the fixed-time parts of planar Wilson loops by weighted sums over 
paths, running from the quark to the antiquark site (a distance R apart). The 
temporal extension of the Wilson loop is denoted as T. In view of the maximal 
spatial distance i? = 6 of quark and antiquark, we have applied Nsmear = 6 
smearing iterations with an optimal smearing parameter (the notation is as in 

(i) 

= 0.125 

= l-4cr""^- (10) 

This parameter had been optimized in order to get an early plateau in T of 
log(W^(i?,T + l)/W{R,T) for the Monte Carlo configurations at /3 = 1.5. For 
simplicity, it has been applied in all measurements of the zero temperature string 
tension. 

The potential is fitted (for all R) from exponential fits over the range T = 
2... 5. For the inversely blocked configurations we have restricted the fit range 
to T = 3. ..5. The potential V{R) has been fitted with a 3-parameter fit that 
includes a constant. Coulomb and linear part over the range R = 1...5. For 
the inversely blocked configurations the Coulomb part is consistent with zero. 
The string tension is obtained from the linear part of the potential. The j3 
dependence and the comparison between normal Monte Carlo configurations 
and inversely blocked (smoothed) configurations is presented in Table y. 

The string tension for inversely blocked configurations has a scaling window 
(assuming asymptotic scaling with the two-loop a(/3)) including our measure- 
ments from 13 = 1.5 to (3 = 1.6. The string tension measured at /3 = 1.8 suffers 
from a small volume effect, since the spatial box size at this (3 is physically not 
larger than 2.35/Tc. For definiteness, we choose the string tension from the 
Monte Carlo configurations in order to express non-perturbatively a(/3) and the 
4-volume of the lattice through the invariant string tension. 

According to the RG philosophy the same large scale physics should be rep- 
resented by the coarse configurations (if obtained by blocking) and the original 
Monte Carlo configurations (on the fine lattice), as long as the blocking scale is 
still safely separated from the physical scale under investigation. It is not clear a 
priori that the inversely blocked configurations U^^ still contain the same large 
scale physics. In case of the string tension this could only be the case if it de- 
couples from the ultraviolet fluctuations. In order to test this it makes sense to 
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compare the string tensions of Monte Carlo and inversely blocked configurations 
at various f3. The string tension is built up at large distances and is expected to 
coincide for both types of configurations, while the perturbative Coulomb part 
of the potential is not expected to survive inverse blocking. We compare in table 
H the zero temperature string tension as found on Monte Carlo and inversely 
blocked configurations, both on the fine lattice. We find our expectation largely 
confirmed. The string tension of inversely blocked configurations is somewhat 
smaller (75 to 90 per cent depending on /3) than the string tension of normal 
Monte Carlo configurations. 

The exception at the lowest /3 = 1.4 can be explained by the fact that at 
very strong coupling the blocking scale interferes with the confinement scale. 
This is supported by the observation that also the topological structure is not 
truly captured. This becomes obvious if one shifts the coarse lattice (in one of 
the 2"* — 1 directions) relative to the fine lattice. This is the reason why we must 
exclude the lowest /3-value also from extracting the temperature dependence of 
the topological susceptibility. 

This is place to compare cooling and inverse blocking as far as the string 
tension is concerned. Cooling conserves the string tension as a function 
of the coohng steps only as long as the topological susceptibility goes through 
a plateau. The topological susceptibility is defined as the plateau value in this 
method . Cooling as used by the MIT group [|T] in order to measure hadronic 
correlation functions goes much farther and docs not conserve the string tension. 

3.2 Is the T ^ String Tension Changing under 
Inverse Blocking ? 

We discuss now our results concerning the temperature dependent string tension 
(t{T). Although this is not our main interest, we want to see to what extent this 
quantity is also insensitive to inverse blocking. The configurations are Monte 
Carlo generated on asymmetric lattices of size 12^ x 4. The string tension (t(T) 
should vanish in the deconfinement phase, at > (3c- A suitable definition of 
the string tension for thermal lattices is provided by the correlator of Polyakov 
lines 



which, in the confinement phase, exhibits the confining potential in the form 



where T^™"-'(x) ~ CT(r)|x| at large distances. V^°^^ (x) contains the self- 
energies of q and q and the perturbative Coulomb potential, too. These are 
expected to change under inverse blocking. 



L(x) ^ tr J|c/xt,4 



(11) 



(L(0)i(x))«exp(-F-7"^(x)/T) 



(12) 
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Figure 2: qq potential from Polyakov line correlators for Monte Carlo (above) 
and inversely blocked configurations at P = 1.5 in the confinement phase 



In the measurement of Polyakov lines and their correlators there is no place 
for smearing and no sense in performing the axial gauge. The normalized cor- 
relators in the Monte Carlo and the inversely blocked ensemble can be directly 
compared for the three /3-values that belong to the confinement phase. Fig. ^ 
shows for f3 — 1.5 the logarithm of the (suitably normalized) Polyakov line cor- 
relator as function of the distance. In Table ^ we collect the extracted values of 
the string tension at the three temperatures in the confinement phase for Monte 
Carlo and inversely blocked configurations. There is an interesting temperature 
dependence of the thermal string tension, but no big difference between the 
string tensions measured on the Monte Carlo and the inversely blocked ensem- 
ble. The reduction is less than 5 per cent at temperatures T < 0.8 Tc and grows 
to more than 20 per cent at T > 0.9 Tc. We observe also at T < Tj, that the 
original Monte Carlo string tension is practically the same as the string tension 
measured on the much smoother inversely blocked configurations and seems to 
be independent of the short range fluctuations. 
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1.40 


1.50 


1.54 


TIT, 


0.656 


0.834 


0.919 


Monte Carlo 
inv. blocked 


0.253(23) 
0.242(8) 


0.193(4) 
0.189(2) 


0.176(5) 
0.137(3) 


Monte Carlo 
inv. blocked 


86(8) 
83(3) 


110(4) 
104(2) 


118(4) 
92(2) 



Table 3: String tension a{T) (above) and a{T)/K\ (below) from Polyakov 
line correlators at various temperatures T < T, in the conEnement phase (with 
two-loop expression for a(f3)AL ) 

4 Gross Topological Features Near to the 
Deconfinement Transition 



The first measurement of topological features of pure SU (2) gauge theory with 
the new method concerns the temperature dependence of the topological sus- 
ceptibility of pure SU (2) gauge theory, here over the range from 0.8 T, to 1.7 T,. 
This question has been studied already in Ref. |^ on lattices of invariable spa- 
tial size (both in lattice units Lg and physical units) by varying the temporal 
extent Lj. We keep the lattice size invariantly equal to 12^ x 4 while varying /3. 
This additional measurement might be of interest as such since now the spatial 
size changes with the temperature. 



4.1 Topological Susceptibility Across the Deconfining 
Phase Transition 

We have measured on Monte Carlo configurations (separated by 20 updates) 
on the lattice 12^ x 4 the topological charges according to Liischer's definition, 
Qinesc/ier [HI |9) , and using the naive topological density The latter is 

defined (actually in a symmetrized way) in terms of plaquettes 

1 +^ 

with 

Qnaive — ^ ^ Qnaive (■^) (-^^) 

X 

(summed over lattice points). The statistics of analyzed Monte Carlo configu- 
rations amounts to 300 configurations. These measurements will be confronted 
for each (3 with corresponding ones on an independent ensemble of 100 inversely 
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blocked (smoothed) configurations. There is no room for details of Liischer's 
charge [|l9|. Important for us is that it is also written as a sum of localized 
contributions 



C luescher 



Qluescher (*^) 



(15) 



(summed over hypercubes). The time consuming part of this algorithm is the 
numerical integration of 2- and 3-dimensional integrals. It should be encouraging 
that evaluating the topological charge of inversely blocked configuration needs 
only of CPU time (function calls) compared to Monte Carlo configurations. 
The topological susceptibility is defined as 



Xtop 



Nsitesa^' 



(16) 



with total charges Q corresponding here always to the naive and Liischer's 
charge. NsUes is the number of lattice points x Lt and we estimate the 



physical lattice volume Ns, 



with the help of the zero temperature string 



tension. Data for Monte Carlo and inversely blocked configurations as function 
of T/Tc (yet without renormalization of the naive topological susceptibility) 



are collected in Table |J. Only the topological susceptibility Xt. 



luescher 
op 



of the 



Xtop 


P 




^luescher / A 4 
Xtop M^L 


^ naive / A 4 
Xtop 1 ^^L 


Monte Carlo 


1.40 


0.656 


0.224(17)10^ 


0.146(12)10^ 




1.50 


0.834 


0.287(24)10'' 


0.370(27)10^ 




1.54 


0.919 


0.310(27)10'' 


0.449(43)10^ 




1.61 


1.089 


0.235(19)10^ 


0.753(59)10^ 




1.80 


1.732 


0.190(16)10^ 


0.308(24)103 


inv. blocked 


1.40 


0.656 


0.127(17)10^ 


0.84(11)10^ 




1.50 


0.834 


0.358(47)10^ 


0.240(29)10^ 




1.54 


0.919 


0.324(43)10^ 


0.239(30)10^ 




1.61 


1.089 


0.138(25)10^ 


0.962(18)10^ 




1.80 


1.732 


0.164(11)10^ 


0.100(70)10^ 



Table 4: Topological susceptibilities in the Monte Carlo and inversely blocked 
samples (with two-loop expression for a(/3)AL J 

inversely blocked ensemble can be interpreted as a topological susceptibility 
(in the sense of the classically perfect topological charge). We have to agree 
that the topological charges belong to the coarse lattice and are defined at a 
scale given by the coarse lattice spacing. The apparent topological susceptibility 
evaluated on the Monte Carlo ensemble of fine lattice configurations is bigger 
than Xtmf''''^^ by one order of magnitude (in the confinement phase) und up 
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to two orders of magnitude (in the deconfinement phase). This does not mean 
that the enormous difference must be attributed to physical fluctuations on the 
scale of the fine lattice spacing. Their contribution can be quantified only if 
one inversely blocks down to the next finer lattice (of size 24'^ x 8) with lattice 
spacing i ■ a. Even lacking this information, the comparison shows how much 
the measurement of the topological susceptibility without the step of inverse 
blocking fails. This means in particular that the (simplified) fixed point action 
alone does not prevent dislocations from contributing to Xtop- 

Let us now discuss to what extent the naive topological charge density op- 
erator can replace Liischer's construction in order to extract the continuum 
topological charge density of a lattice configuration. It is known for simula- 
tions with Wilson's action in the range f3 — 2. to 3., that the corresponding 
perturbative renormalization factor |32[ is very small compared to one. Also 
for simulations with the fixed point action, the topological susceptibility, eval- 
uated immediately for Monte Carlo configurations in the confinement phase, is 
two orders of magnitude smaller with the naive charge than the susceptibility 
defined by means of Liischer's charge Qiuescher ■ 



a 




i.O -6.0 -4.0 -2.0 0.0 2.0 
Q_naive 



6.0 8.0 10.0 



Figure 3: Scatter plot of Liischer's vs. naive topological charge for inversely 
blocked configurations at jS = 1.5 

In contrast to this, the susceptibihties Xtop^'^ ^^'^ xiZp^'^'^'^^ differ by less 
than 30 per cent for the inversely blocked ensemble of configurations. The 
remaining discrepancy can be illustrated with Fig. 0. This scatter plot of 
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inversely blocked configurations shows how Qiuescher is correlated with Q naive- 
There is only a small variance between the two charges evaluated after one 
step of inverse blocking. The naive charge differs from Liischer's charge never 
by more than one unit. We can define a renormalization factor for the naive 
charge Zq^\p) on inversely blocked configurations by the slope of the scatter 
plot. The corresponding numbers are collected in Table ^ for our set of /3- values 
near to the phase transition.^ The renormalization factor can be used to relate 
the continuum to the naive topological density after the first step of inverse 
blocking: 

^(^)w(^)a' = gile(^)«' • (17) 





1.50 


1.54 


1.61 


1.8 




.826 


.877 


.843 





Table 5: Effective renormalization Zq factors of the naive topological charge 
density for one step of inverse blocking at various values of j3 on the 12^^ x 4 
iattice 

Further steps of inverse blocking will bring zjj'^\l3) nearer to one. Com- 
pared to the renormalization factor for Monte Carlo configurations, this effect 
is similar to the effect of improving the operator of topological density |^ pro- 
posed by Di Giacomo et al. These operators are constructed exactly as in ( p^ 
in terms of suitably smeared links C/*^"-' (iteratively smeared to the n-th level). 
This or other improved operators of topological density |3^, ^ can be used 
for analysing inversely blocked configurations. For the topological susceptibility 
the renormalization factor as derived from the scatter plot is sufficient. When 
the naive charges (that went into the susceptibility of smoothed configurations 
in Table are multiplicatively corrected by the Zq^\l3) factors from Table || 
we get a unique /^-dependence of Xtop- We have shown in Fig. ^ the topolog- 
ical suceptibility for our four temperatures. The lattice spacing a(/3) has been 
non-perturbatively expressed through the zero temperature string tension as 
measured at each respective /3-value on the 12^ lattice. 

Assuming the value a = (440MeV)^ for the zero temperature string tension 
we estimate at T = 0.834 Tc a susceptibility Xtop = (165.5MeV)l Compared 
with the standard value of Xtop — (180MeV)^ at zero temperature there is not 
much room for increase at smaller temperature. 

For building instanton models of vacuum structure it is important to know 
not only the topological susceptibility but also the average density of instantons 

^For our highest /3 non-zero charges are very rare, such that the Z^^^ factor is not reliably 
known. 
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Figure 4: Topological susceptibility for inversely blocked conEgurations across 
the deconGning phase transition for Liischer's topological density 

and antiinstantons. Only in the most simple-minded dilute gas model the topo- 
logical susceptibility coincides with this density. In addition to the total charge 
(eqs. or (|l^)), the following operator is tentatively chosen to represent 

7V+ + , the number of instantons plus antiinstantons 



We can understand this, somewhat loosely, as a measure for the temperature 
dependent average glueball field (gluon condensate). Let us call it "topologi- 
cal activity". Compared with measurements of this quantity on Monte Carlo 
configurations, it is reduced by ^ to ^ on inversely blocked configurations, 
depending on temperature. For normal Monte Carlo configurations there is too 
much noise at low values of Without inverse blocking, Atop is unreli- 

able as an estimator of the number of instantons plus antiinstantons. In Fig. 
^ we show the average topological activity of inversely blocked configurations 
as function of T/Tc- Between (3 = 1.54 and (3 — 1.61 (over temperatures from 
T = 0.92 Tc to T = 1.09 Tc) it drops by 50 per cent, but the decrease becomes 
slower at higher temperature. 

*In contrast to that, the topological activity measured on Monte Carlo configurations is 
less suppressed as a function of temperature. 
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Figure 5: Average topological activity (Atop) (according to Liischer's charge) of 
inversely blocked configurations as function ofT/Tc 



Concentrating now on the inversely blocked configurations, the topologi- 
cal susceptibility is strongly suppressed with increasing temperature, compared 
with the decrease of the topological activity. Fig. |6| demonstrates this additional 
suppression of uncompensated topological charge in the deconfinement phase. 
Apart from this observation, an interesting low temperature bound is suggested 
(see the dotted horizontal line in Fig. ||) : 

6 (Atop) 6 

There is a factor ^ = jjjf- (equal to the entropic bound) in this relation where 
b is the one-loop coefficient in the QCD /3-function. Models describing this 
0{1/Nc) suppression of the topological susceptibility compared to the densities 
n± of instantons and antiinstantons have been discussed in an instanton liquid 



model with scale invariant hard core interaction |36 and in a similar droplet 
model Hi. 

We have found that (i) the topological susceptibility for smoothed configu- 
rations is of reasonable size on the confinement side of the transition, (ii) the 
topological susceptibility is smaller by a factor | than the topological activity 
at low temperature and that (iii) the topological susceptibility decreases more 
than the topological activity entering the deconfinement phase. 
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Figure 6: Ratio of topological susceptibility xtop to average topological activity 
{Atop) / N sitesCi^ (according to Liischer's charge) of inversely configurations as 
function ofT/T^ 



4.2 Monopole Content of Monte Carlo and Inversely 
Blocked Configurations 

The Abelian monopole degrees of freedom are exposed by putting the lattice con- 
figurations into the maximally Abelian gauge . Only few investigations have 
addressed the question how the dynamics of Abelian monopoles gets modified 
at the deconfinement transition and within the deconfined phase Abelian 
dominance is expected to become manifest at large distances. Hence it is nat- 
ural to study these aspects in the ensemble of inversely blocked configurations 
which are smoothed at small distances. 



It is not necessary to recall here details on the Abelian projection |13|. By 
a gauge cooling algorithm configurations are iteratively put into this gauge, 
whether they are directly Monte Carlo generated or obtained by inverse block- 
ing. Finally, each configuration is mapped to a corresponding Abelian one. 
The links are represented by U{1) phase factors and the non-diagonal fields are 
dropped. Using this new ensemble of configurations, confinement can be studied 
either in terms of Abelian Wilson loops formed out of the Abelianized links or, 
alternatively, in terms of the monopole currents. These are detected JTi] seeking 
for magnetic flux penetrating the surface of a 3-dimensional cube. Each time a 
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non-zero flux is detected, the dual (with respect to the cube) link x,^ is said 
to carry a monopole current m^,{x). We do not discuss here the correct descrip- 
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Figure 7: Ratio between the averages of total monopole current number M 
and topological activity Atop as function ofT/Tc for Monte Carlo and inversely- 
blocked configurations 



tion of monopole condensation |3S 
the monopole current numbers 



\'m^,[x)\ 



and restrict ourselves to measurements of 



(20) 



It is important to distinguish spacelike (/i = 1,2,3) and timelike (/i — 4) dual 
links carrying monopole currents. We will see that their numbers behave differ- 
ently at the deconfinement temperature. 



There are indications that the monopole current number = X^u -^^m 



IS 

related to the topological activity Ato-p- Both quantities are reduced by the 
same factor ^ to at the lowest and the highest temperatures, respectively, 
as the result of blocking followed by inverse blocking. We show in Fig. ^ 
the ratio {M)l{Atop)- This ratio is a smoothly falling function of temperature 
across the phase transition for normal Monte Carlo configurations. For inversely 
blocked configurations, however, we find an sudden drop of this ratio at Tc- This 
is mainly due to the suppression of spatially directed monopole currents. The 
corresponding ratio for timclikc monopole currents alone (see Fig. ^ is a smooth 
function of temperature for Monte Carlo and inversely blocked configurations. 
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Figure 8: Ratio between the averages of timelike monopole current number Mt 
and topological activity Atop as function ofT/Tc for Monte Carlo and inversely 
blocked configurations 

The ratio between Mg (summed over = 1,2,3) and 3 Mt (/i = 4) is ex- 
pected to change at the deconfincmcnt temperature In Fig. ^ the ratio 
(Ms)/3(Mf) is shown for Monte Carlo and inversely blocked configurations as 
function of T/T^. For the original Monte Carlo configurations the monopole 
current number is space-time symmetric in the confinement phase and the ra- 
tio (Ms) /3{Mt) changes only slightly over the temperature range interpolating 
from the confinement to the deconfining phase. The deconfinement effect is not 
strong because it is hidden by the huge overall monopole activity in the Monte 
Carlo configurations. The effect becomes clear only for inversely blocked con- 
figurations. There exists an surplus of timelike monopole currents already in 
the confinement phase (which probably must be interpreted as a finite Lt effect 
produced in the process of blocking). There is a strong reduction of space- 
like monopoles at T = 1.089Tc, and deeper in the deconfined phase they are 
completely suppressed. 

Fig. |l^ expresses the monopole currents for inversely blocked configurations 
in the form of 3-dimensional densities 

^ {M^) (21) 



Nsitesa^ 

of timelike and spacelike monopoles (in units) as a function of T/T^- Mo- 
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Figure 9: Ratio between the averages of timelike and spacelike monopole cur- 
rents across the deconhning phase transition for Monte Carlo and inversely 
blocked configurations 



nopoles, condensed in the confining phase, become static (massive) particle- 
Hke objects in the deconfined phase with a density rising with temperature 
(cf. Ref. ). One may speculate about their relation to 't Hooft-Polyakov 
monopoles. Spacelike monopole currents (responsible for confinement below 
Tc) become strongly suppressed at Tc and have completely disappeared at T ~ 
1.7 T,. 



5 Resolving the Vacuum Structure 
5.1 Topological Density Correlation 

Information on the nature of topological excitations is contained in the point- 
to-point correlation function of the topological charge density. This is easy to 
measure, but for usual Monte Carlo configurations [7*^"^ this measurement has 
no value. 

Several groups have attempted to define shape and size of topological ex- 
citations by the use of cooling techniques. The MIT group ^ has used for 
this purpose deeply cooled configurations. There is an uncertainty about the 
loss of essential structures during cooling. Sometimes there appears a partic- 
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Figure 10: S-dimensional densities of timelike and spacelike monopoles across 
the deconGning phase transition for inversely blocked conBgurations 



ular structure just during the first few cooling steps |16, M. Recently, there 



have been attempts to improve cooling by choosing suitably improved actions 



[ p3| , 34, |35[] for (unconstrained) relaxation. Partial success has been achieved 
to stabilize instantons above some threshold in size. But certain topological 
structures (instanton-antiinstanton pairs) unavoidably disappear in the process 
of cooling. The method of inverse blocking is welcome because it enables to 
detect correlations on a set of well-defined configurations, which are smooth 
but interpolate, in a locally correct way, coarser equilibrium configurations. 

In Fig. |ri| we compare the normalized correlation function of the topological 
charge density measured on the 12^ x 4 lattice in the confinement phase at 
P = 1.5 on Monte Carlo and inversely blocked configurations. For normal Monte 
Carlo configurations there is no signal at any distance except for x = 0. In order 
to demonstrate the influence of temperature on the signal we show in this figure 
also the correlation function measured at /3 = 1.8 in the deconfinement phase. 
Due to the fixed normalization at a; = only the change of the instanton radius 
can be seen (as explained below). 

Smoothing is necessary in order to expose a signal in the density-density 
correlation function not only in the case of the naive topological density def- 
inition qnaivei^)- Fig. |l^ shows the same for Liischer's charge density. The 
effect is similar, but the correlation of the naive charge density always extends 
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Figure 1 1 : Normalized correlation function of the naive topological charge den- 
sity at /3 = 1.5 (for Monte Carlo and inversely blocked configurations) and at 
(3 = 1.8 (only inversely blocked). Solid lines are fits according to eq. (22), (23) 



somewhat farther. This can be understood due to the extended nature of the 
operator (|l3|). 

Getting a reasonable fit of the correlation function to the analytical shape 
derived from the instanton profile 

q^nstix) = (^^) (22) 

by folding 

{q(x)q(Q)) oc qinst jx - z)qinstiz), (23) 

z 

we have determined the average instanton radii pinst as given in Table |^. The 
corresponding fitting curves are shown in Figs. ^ and|l^. The instanton radius 
as determined by this method decreases more rapidly with rising temperature 
than the inverse temperature. Some change of Pinst at deconfincment has been 
observed by the cooling method, too |^. 
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Figure 12: Normalized correlation function of the geometric Liischer topological 
charge density (as in Fig. 11) 

5.2 Monopole— Instanton Correlation 

It has been demonstrated before for certain classical configurations that there 
is a close connection between instantons and monopole world loops. One can 
visualize this in the maximally Abelian gauge |^5j . A nontrivial correlation has 
been measured between topological charge density and monopole currents in 
the Monte Carlo ensemble of Euclidean field configurations, too We show 
here in Fig. the analogous normalized correlator for the original Monte Carlo 
ensemble (generated with the fixed point action at f3 — 1.5 in the confinement 
phase) and for the corresponding ensemble of inversely blocked configurations. 
The latter correlator (for smoothed configurations) is slightly wider. In contrast 
to the topological density-density correlations, this correlator can be detected 
in the Monte Carlo configurations as well. 

Comparing this with the same (normalized) correlation function at P — 1.8 
in the deconfined phase (not shown here), we find that they are identical as a 
function of the distance in lattice units within the error bars. This means that 
this correlation length changes proportional to the inverse temperature. 
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p 


1.50 


1.54 


1.61 


1.80 


Liischer's charge 
naive charge 


1.4 a 
1.8 a 


1.4 a 
1.8 a 


1.35 a 
1.75 a 


1.1 a 
1.6 a 



Table 6: Instanton radii pinst in lattice units as estimated from the topological 
charge density correlation function at two (3 values in the confinement and two 
f3-values in the deconfinement phase, according to the two charge definitions 
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Figure 13: Normalized correlation function of monopole currents with topolog- 
ical charge density (according to Liischer) at (3 = 1.5 in the confinement phase 

6 Conclusion and Outlook 

We have explained in this paper how inverse blocking can be used to analyze 
topological structures of genuine quantum configurations, usually hidden under 
ultraviolet fluctuations. 

Apart from cooling techniques, inverse blocking is the only way to measure 
autocorrelations of the topological density giving information on the size of in- 
stantons. Information on monopole currents is available also for normal Monte 
Carlo configurations. Monopoles are correlated with the topological density, 
both in Monte Carlo and inversely blocked configurations. However, the den- 
sity of monopole currents is reduced, roughly proportional to the topological 
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t=4/MQ_V0R.DAT 




Figure 14: Liischer's topological charge density for \q(x)\ > go aiid monopole 
world lines in a timeslice of a typical Monte Carlo configuration in the confine- 
ment phase; signs of q are distinguished by different color 

activity, if a configuration undergoes the procedure of blocking followed by in- 
verse blocking. 

Only a small part of Abelian monopoles seen usually in Monte Carlo con- 
figurations seems to be really important for confinement. This part survives 
smoothing in the confinement phase. We come to this conclusion because the 
string tension is found to be largely unaffected by the removal of ultraviolet fluc- 
tuations. This observation supports semiclassical scenarios of confinement. In 
the deconfinement phase, however, we find that spacelike monopoles are strongly 
suppressed for inversely blocked configurations while only timelike ones survive. 

Armed with inverse blocking, we can do better than measuring the topo- 
logical density-density correlation or the correlation between monopoles and 
topological density. We are now able to study clustering properties of the topo- 
logical charge using inversely blocked lattice configurations. We can relate clus- 
ters of topological charge to the diluted network of monopole currents. Instead 
of examining idealized semiclassical configurations under different aspects of 
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topological and monopole structure, the method of inverse blocking makes it 
possible to analyze generic Euclidean lattice field histories in order to charac- 
terize the physical phase one is dealing with. The results should be interesting 
for vacuum model builders. In view of chiral symmetry breaking and restora- 
tion, it would be also interesting to study the influence of the inverse blocking 
RG transformation on the spectrum of the Dirac operator. 

The study of clusters of topological charge promises to give a more direct 
access to the parameters and correlations that should be built in an instanton 
model for QCD at zero and finite temperature. In Figs. |lj and ^ lattice 
hypercubes of one timeslice are marked if the Liischer topological charge density 
exceeds a threshold value qo = 0.015. Different colors symbolize the sign. The 

t=4/MQ_NACH.DAT 




Figure 15: Clustering of the Liischer topological charge density and monopole 
world lines in a timeslice of the same configuration as in Fig. ^ but after 
smoothing 

original Monte Carlo configuration does not allow to guess any structure. The 
other picture shows the structure emerging after blocking and inverse blocking. 
The effect is the same using either naive or Liischer's charge density. Therefore 
one can abandon to use complicated algorithms and use naive or improved loop- 
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oriented definitions of topological charge density. 
There are two immediate directions for further studies. 

• Study of the clustering properties of topological charge in individual lattice 
configurations and of the relation to the monopole currents. 

• Repeating the present explorative study at higher (3 with a particular 
emphasis on the scale dependence of topological structure. 

Further progress in the construction of the perfect action should be incorporated 
into this program. 

Note added: 

During the time elapsed since the first submission of our article, several in- 
teresting papers have appeared which have attempted to investigate the space 
time topological structure of lattice gauge fields. Three methods have been 
used: (improved) cooling, fermionic (spectral) methods and smoothing. The 
latter method had been proposed in this paper as a general, renormalization 
group oriented pattern recognition method for large scale structures in terms of 
topological charge and monopoles which conserves the string tension. 

While the subsequent studies were mostly focussed on vacuum (zero tem- 
perature) configurations, the results are still controversial as emphasized in Ref. 
[ p3| . As a generic rule, one can state that improved cooling is biased towards 
larger instantons (because only those are stable) while the topological suscep- 
tibility tends to be lower; in contrast, the smoothing method (as developed 
further by the Boulder group |40| ) favours smaller instanton sizes and tends 
to find somewhat higher topological susceptibility. This applies to SU (2) pure 
gauge theory where the results can be compared (but not without additional 
assumptions). 

We want to stress an important difference of Ref. |Q to the work published 
here. Smoothing has been applied there in an iterative way that leads finally to 
locally classical configurations. The number of iterations appears there as some- 
what subjective parameter similar to the number of (unimproved or improved) 
cooling steps. 

This goes far beyond our intention. From our point of view, the advantage 
of smoothing is that it provides a minimal smoothing of quantum fluctuations in 
clear correspondence to a certain scale. This scale is well-defined by the highest 
blocking level. 

The effect of smoothing on various Abelian projections has been critically 
examined in Ref. The maximally Abelian gauge is the only one which 

seems to exhibit long range physics also if applied to smoothed lattice fields. 
While we had concentrated on the density of monopoles, it has been found there 
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that the property of Abelian dominance (i.e. dominance of the Abehan string 
tension) is only partly preserved. 

A large number of investigations was studying topological structures using 
fermionic methods going back to Smit and Vink p2) . It is important in the 
context of this work that these methods seem to work immediately, without 
smoothing or cooling techniques which otherwise would have to be applied to 
the Monte Carlo lattice gauge field configurations. Investigating the spectral 
flow of the lowest eigenvalues of the Dirac operator |4j, ^ these authors 
are able to count the number of localized lumps (and identify the respective 
sign) of topological charge. The local topological structure is thereby indirectly 
accessible through the localization properties of the corresponding eigenmodes. 
This has been checked in Refs. ^ with cooling, but there is no comparison 
so far with the topological structure detected by smoothing. 

We have recently improved our relaxation algorithm for implementing inverse 
blocking [1^. We are now in the position to optimize a fixed point action 
of moderate complexity for a given inverse blocking parameter k. We could 
check a new parametrization of the perfect action for SU{2) which is truncated 
compared with the action used in Ref. |Q. We thank the Boulder group for 
correspondence about that and for providing the truncated action [Q . We can 
now show that the parametrization of Table |l| is not compatible with the k — 12 
required. This inconsistency had forced us to accept (in this paper) a smaller 
effective k. Wc have convinced ourselves that the statistical observations are 
independent of this choice although the inversely blocked configurations may 
differ locally. 
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